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Some  integrals  are  presented  that  can  be  expressed  in  terms  of  the  QM-ftmction, 
which  is  defined  as  /  2  2\ 

m*/.)"'1  «p(-  . 

and  where  Ijg-i  is  tee  modified  Bessel  function  of  order  M-l.  Also,  some  integrals 
of  tee  C  jf -function  are  evaluated.  All  derivations  are  included  in  tee  appendix. 
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SOME  INTEGRALS  INVOLVING  THE  QM-FUNCTION 


INTRODUCTION 

J 

The  detection  probability  for  a  nonfading  signal  present  on  M  multiple 
observations  and  subjected  to  £  quadrat  e  combination  of  recover  matched- 
filter  outputs  is  given  by  the  Qj^-function^Ref.  l] : 


«  M-l  /  2  2\ 

QM(«.l»S  exp(-i-p-)lM1(ax), 


M>  1 . 


This  is  also  called  the  noncentral  chi-square  distribution  [Ref.  2] .  Physically, 
Qm  is  the  probability  that  the  random  variable  k|t  [<*k  +  °k)2  +  &k  +  dk>2^ 
exceeds  threshold  T ,  where  {xfc|  and  are  zero-mean  independent 
Gaussian  random  variables  with  Identical  variances,  *2,  and  {c,K}  and 
{d^}  are  arbitrary  constants.  The  parameters  a  -  [  £(  (c2  +  d^)/*2]1'2 
and  b=T/v.  \ 

When  the  signal  amplitude  is  subject  jto  hiding,  (1)  must  be  averaged  with 
respect  to  the  distribution  .of  such  fadings.  In  ids  report,  several  integrals  of 
the  Qj^-function  wtH-be^  evaluated  that  aid  in  computing  the  performance  in 
fading  channels;  direct  numerical  integration  of  die  QM-function,  which  is 
time-consuming,  is  therefore  circumvented.  Numerical  evaluation  of  QM  is 
considered  in  Re&<£.  jl  ^  •  \  v  , 


A 


GENERAL  RELATIONS  FOR  QM-FUNCTION 
Derivations  of  the  succeeding  relations  are  presented  in  the  appendix: 


V°’b>  ■  ••(- t)  E  ff  (t)  •  Q(-b» ' 1  • 
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/bvM-1  /  2  +  fa2\ 

V‘’b)  -  Q^o.b)  +(;)  «p^-  —5—  J  iM.i<llb) 

/  2  U2\M^  ...k 

- Q(a,b)  +  «p(- (l)  Ik<a^*  Mi2  '  ^ 

'  'k«i 

where  Q  =  Qi  ie  the  special  case  M  =  1.  The  relation  Qo(a,b)  =  1  -  Q(b,a) 
also  follows  from  the  first  line  of  (3)  and  Ref.  3,  (5).  In  the  following  we  will 
generally  consider  the  integer  M  >  1 ;  however  many,  but  not  all,  of  the  fol¬ 
lowing  equations  apply  also  to  M  =  0 . 

QM<a,b)  =  exp(-  £2  (|)  yab)  .  <4 


k— M+l 


QM<a’b)  =  1  -  «P 


(•=¥)£©“** 

k=M 


'  •  k--M+l 


QM(a,b>+QM<b,a)  =  1  +  C*P 


'  k--M+l 
M-1  . 


-fj1  *  +  2  2  ’k^]!  • 

1  k-1 


j(a,b)  =  exp^- 


2\  /  2/  /  2/  y* 

k«0  n«0 


aQ„(a,b) 


dQM(a.b) 

.  j 

(,b\M'1  t 

r  2  .2\ 

'  a  +b 

a* 

i :  exp' 

f  2 
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The  numerical  evaluation  of  Q^(a,b)  considered  in  Ref.  2  is  based 
upon  a  modification  of  (8).  Notice,  from  (8),  that  QM(a.b)  is  an  even  function 
of  both  a  and  b  for  a  and  b  real.  Also,  from  (1),  Qj^(a,b)  is  an  analytic 
function  of  a  and  b  for  all  finite  complex  a  and  b. 


INTEGRALS 

Numerous  other  integrals,  in  addition  to  those  listed  below,  are  available 
by  the  use  of  (3)  and  Ref.  3,  provided  the  integrals  on  the  Bessel  functions  in  (3) 
can  be  evaluated  in  closed  form.  Rather  than  exhaustively  list  all  these  possibili¬ 
ties,  some  of  the  simpler  integrals  are  considered  and  evaluated.  Extensions 
to  other  integrals,  which  may  be  of  interest  to  die  reader,  may  profitably  be 
pursued  via  (3),  Ref.  3,  and  a  good  table  of  integrals  such  as  found  in  Ref.  4. 

In  the  integrals  below,  the  parameters  are  presumed  to  be  real  and  posi¬ 
tive.  The  results  may  be  generalized  to  negative  or  complex  parameter  values 
in  many  cases  by  symmetry  or  analytic  continuation.  The  integrals  have  been 
checked  numerically. 

An  integral  closely  related  to  QM  is  defined  by 

/  g  2\ 

dx x(a)  (U) 

b 

Some  of  the  properties  of  this  function  follow: 

W‘'h|  *  qM(a'b)  -  (!)  “p(-  5-rL)IM(ab>  •  (12) 


qQ(a.b)  =  Q(a,b)  . 


/  2  2\  M-l  .jj 

,(a,b)  -  Q(a,b)  -  exp(-  S-ti-j  £  (ij  yah)  . 


qM<a,b)  *  1  '  QM*b,a)  ’ 


lM(l1' *>  *  2  I1  '  “K  W  ["o^2)  +  2  2  ! 

1  L  k-1  J7 


(13) 

(14) 

(15) 

(16) 
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BESSEL  FUNCTIONS,  EXPONENTIALS,  AND  POWERS 

Jdx  exp(-p2x2/2)  JM_j(«)  JM(bx)  =  -(g)  [l  -  QM<a/P*  b/P>]  • 
0 

^dx  exp(-p2 x/2)  JM_1(ax) 


(17) 


(18) 


in  integral  (20)-(21).  .=  'JiT#  .  »  -  V«(P  -  «>  .  *  *  V«(P  ♦  •>  • 


^dx  exp(-px) p>b'  <20> 

L  «xp(-p=0  -  .)M[l  -«M(U.  ’O] -(P+.)M[l- «„(».»)]}•  P'1"- 

0  ,b  (21) 

{*  *“  «p(-pV/2) :„.!<«)  *^M  «M<a/Pl  bp)  • 

b 

^  X-M+1  «xp(-p2x2A)  !„(“)  -  “P^)  [1-«M<bP’a/P)]  • 


4 


(23) 
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Q-FUNCTION,  EXPONENTIALS,  AND  POWERS 


where 


k  <  N  -  1 


1  +  ^  ,  k  =  N  -  1 

a 


Qm-FUNCTION  and  powers 


ID 

(dxx* 


1+M 

2 


°2  r(M+i^.) 


Q«(b.ax)  - 
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Q  -FUNCTION,  EXPONENTIALS,  AND  POWERS 
M 

J  dx  x2M_1  exp(-p2x2/2)  QM(ax,b) 


(dxxp  expf-pV/z)  Q  (ax,b)  =  G 
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dx  x  exp(-p2x2/2)  QM(ax,b)  =  —  exp(-p2c2/2)  QM(ac,b) 

P 
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Appendix 

DERIVATIONS 


(2)*:  The  first  two  relations  in  (2)  follow  directly  from  definition  (1); 
Ref.  4,  6. 631  4;  and  Ref.  5,  9. 6. 3.  The  third  relation  results  if  we  employ 
Ref.  5,  9.6.7,  and  the  integral 


*°  M-l  -x  Mil  k 
x  e  -c 


(M  -  1)1 


(A-l) 


k-0 


The  fourth  relation  follows  by  use  of 

b 


va'b) = 1  -  5 dx  x(f)  “"l-  H1-) 


(A-2) 


and  as  a  result  of  the  asymptotic  behavior  of  IM1  (Ref.  5,  9. 7. 1). 

(3) :  Integrate  by  parts  with 

/xXM-l  /  2  2\ 

u  =  (a)  =  (A"3) 

and  employ  Ref.  5,  9. 6. 28.  The  second  line  in  (3)  follows  by  repeated  applica¬ 
tion  of  the  first  line. 

(4) :  Use  (3);  Ref.  3t  (3);  and  Ref.  5t  9. 6. 6. 

(5) :  Use  (3)  and  Ref.  3,  (4;. 

(6) :  Use  (3);  Ref.  3,  (5);  and  Ref.  5,  S.  6. 6. 

(7) :  Let  b  =  a  In  (6). 


♦The  number  at  the  beginning  of  each  paragraph  refers  to  the  equation  of 
that  same  number  in  the  main  text  of  this  report 
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(8) :  Expand  IM-1  in  (1)  in  a  power  series  and  integrate  term  by  term 
ustig  (A-l). 

(9) :  Differentiate  (1)  with  respect  to  b . 

(10) :  Differentiate  (?)  with  respect  to  a,  and  employ  Ref .  3,  (8),  and 
Ref.  5,  9.6.28. 

(11) :  This  is  a  definition  of  qj^(a,b). 

(12) :  Integrate  (11]  by  parts  using 

/  \"M  /  2  x  2\ 

-(;)  hP*-  <A'4) 

and,  from  Ref.  5,  9. 6. 28, 

£{(£)  y«y)}-Ki)  W“y)  •  <A-5> 

Then  use  (11)  again. 

(13) :  Use  definition  of  Q(a,b)  in  (1). 

(14) :  Use  (12)  and  (13). 

(15) :  Add  (14)  to  the  second  line  of  (3),  with  a  and  b  interchanged; 
use  Ref.  3,  (5). 

(16) :  8et  b  *  a  in  (15)  and  use  (7). 

(17) :  From  Ref.  4,  6.633  2, 

0 

<A'6) 


a  a 
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the  last  step  was  obtained  via  Ret  5,  9. 1. 30.  Integrate  (A-6)  with  respect  to 
b  and  use  (1)  and  (2). 

(18) :  Set  b  =  a  in  (17)  and  use  (7). 

(19) :  Replace  a  by  1a  in  (18)  and  use  Ref.  5,  9.6.3  and  9,6.30. 

(20) :  Utilize  the  method  used  for  (24)  in  Ref.  3,  p.  24;  employ  (4)  and 
(5)  to  sum  the  series;  and  then  use  (3). 

(21) ;  Use 

Jdx  g(x)  =  jdxg(x)  dx  g(x)  (A-7) 

0  0  a 


conjunction  with  (2)  and  (20).  The  special  case  of  p  =  b  is  available  directly 
Ref.  5.  11.3.13. 

(22) :  Use  definition  (1). 

(23) :  Employ  (5)  and  (15). 

(24) :  Express  Q  in  integral  form  via  (1)  (for  M  =  1) ;  Interchange  inte¬ 
grals;  and  use  (A-l),  Ref.  4,  6.631  1,  and  Ref.  5,  13.1.27  and  13.6.9.  1^ 
is  a  Laguerre  polynomial. 

(25) :  Use  Ref.  3,  (5);  (A-l);  Ref.  4,  6.631  1;  and  Ref.  5,  13.1.27  and 
13. 6.9. 

(26) :  Use  (1);  Interchange  integrals;  and  use  Ref.  4,  6. 631  1. 

(27) :  Use  (26)  and  Ref.  5,  13. 6.9. 

(28) :  Express  QM  in  integral  form  via  (1);  interchange  integrals;  and 
use  Ret  4,  6.631  4,  and  (A-l). 
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(29) :  Employ  the  upper  line  of  (3)  and  Ref.  4,  6. 631  1.  For  p  an  odd 
integer,  the  recursion,  (29),  can  be  started  via  (Ref.  5,  chap.  13): 

^(OjMjx)  =  1 , 

I"  M-2  k 

^(ljMjx)  -  (M  -  1)!  x  exp(i)  -  £  ~  , 

L  k=0  . 

1F1(n;MpO=^~|(M+l-n)  1F1(n-2;Mpc)  +  (2n-2-M+x)  ^(n - l:Mpc)J ,  n>2  . 

(A-8) 

For  p  an  even  integer,  we  need  iFi(f  “d  recursion. 

(30) :  Use  (1);  interchange  Integrals;  employ  Ref.  4,  6,631  1,  and  (A-8) 

and  (A-l);  interchange  the  double  sum;  and  employ 


n-1 


XX-V^r- 


(A-9) 


k-0 


(31):  Express  QM  in  integral  form  via  (1),  interchange  integrals,  and 
use  (22)  twice.  If  we  replace  p  by  ip  in  (31),  the  integral  converges  for  p  <  a. 


(32):  Integrate  by  parts,  with 

u  *  QM(ax,  b),  dv  =  dx  x  exp  (~p2x2/ 2)  , 


(A-10) 


and  use  (10)  and  (23). 

(33):  Take  the  partial  derivative  of  the  integral  with  respect  to  b ;  use 
(9)  and  Ref.  4,  6. 633  2;  then,  utilising  the  fact  that  the  integral  in  (33)  is  zero 
at  b  *  «o  ,  integrate  on  b  and  use  (22). 


TR  4755 


REFERENCES 

1.  C.  W.  Helstrom,  Statistical  Theory  of  Signal  Detection,  second  edition, 
Pergamon  Press,  New  York,  1968,  section  VII,  eq.  (2. 18). 

2.  G.  H.  Robertson,  "Computation  of  the  Noncentral  Chi-Square  Distribution,  ’’ 
Bell  System  Technical  Journal,  vol.  48,  January  1969,  pp.  201-207. 

3.  A.  H.  Nuttall,  Some  Integrals  Involving  the  Q -Function.  NUSC  Technical 
Report  4297,  17  April  1972. 

4.  I.  S.  Gradshteyn  and  I.  W.  Rhyzik,  Table  of  Integrals,  Jertes,  and  Products, 
Academic  Press,  New  York,  1965. 

5.  Handbook  of  Mathematical  Functions.  U.  S.  Department  of  Commerce, 
National  Bureau  of  Standards,  Applied  Mathematics  Series  No.  55,  June 
1964. 


13/14 

REVERSE  BLANK 


